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A measurement of the ti production cross section in pp collisions at ^/s =1.96 TeV using events 
with two leptons, missing transverse energy, and jets is reported. The data were collected with the 
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CDF II Detector. The result in a data sample corresponding to an integrated luminosity 2.8 fb 
is: 

fftt = 6.27 ± 0.73(stat) ± 0.63(syst) ± 0.39(lum) pb. 
for an assumed top mass of 175 GeV/c 2 . 



PACS numbers: 14.65.Ha, 12.38.Qk, 13.85. Qk 

I. INTRODUCTION 

This paper describes a measurement of the it 
production cross section in pp collisions at y/s 
= 1.96 TeV with the CDF detector at the Fer- 
milab Tevatron. This measurement requires the 
identification of both leptons in the decay chain 
it -> (W+b)(W-b) -> (£+v e b)(e~v e b). Events 
are selected with two high transverse energy lep- 
tons, high missing transverse energy (Ifir) and at 
least two jets in the final state. From the excess 
of events selected in the data over the predicted 
background from other known standard model 
(SM) sources we obtain a measurement of the 
production of it events. 

The top quark pair production in the standard 
model proceeds primarily by quark-antiquark an- 
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nihilations. At the Tevatron the predictions 
are 85% quark-antiquark annihilations and 15% 
gluon fusions. At the Large Hadron Collider at 
y/s = 14 TeV the situation is predicted to be 
very different with 90% of the production being 
due to gluon-gluon fusion and 10% due to quark- 
antiquark annihilation. 

This analysis improves upon a previous mea- 
surement of the cross section using the same 
dilepton (DIL) selection in a data sample with 
an integrated luminosity of 0.197 fb -1 [l|. Unlike 
other CDF measurements of the it cross section 
in the dilepton channel 0], where one i is iden- 
tified as 6 or while the other is identified by 
the presence of a high momentum central track, 
the DIL analysis positively identifies both leptons 
as either electrons or muons from W decays or as 
products of semileptonic decays of r leptons, thus 
allowing for the comparison of the observed yield 
of it decays to ee, fj.fi and e/i final states with the 
predictions from lepton universality. 

The measurement provides a test of the QCD 
calculations of the it cross section Q in a chan- 
nel which is independent and complementary to 
other measurements of the it cross section in 
higher statistics final states where at least one 
W boson from the top quark is reconstructed 
via its hadronic decay, W — > qq'. The dilep- 
ton final state suffers from a lower statistical 
precision, as the product of the branching ra- 
tios of the semileptonic W decay BR(W + — >• 
£+v)xBt(W- -> i~v) pb 5% with t = e or (i, 
but it has a signal to background ratio well above 
unity even before requiring the identification of 
one of the jets originating from a b quark. This 
analysis does not require jets in the events to have 
secondary vertexes consistent with the presence 
of a 6-hadron decay as this selection would fur- 
ther reduce the acceptance by almost 50%. 

In Sec. |TT] we give a short description of the 
detector. In Sec. IIIII the data sample and event 
selection are presented. Section ITVl presents the 
formula used for the cross section calculation and 
the measurement of the it acceptance in dilepton 
events. Section [V] describes the calculation of 
the backgrounds. Systematics uncertainties are 
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covered in Sec. IVII In Sec. IVIII observations are 
compared to predictions in control samples char- 
acterized by the presence of two leptons plus high 
IfiT in the final state. We conclude by presenting 
the result of our measurement in Sec. lIXl followed 
by a short summary in Sec. |Xj 



ing and thus determine the beam luminosity The 
total uncertainty on the luminosity is estimated 
to be 5.9%, of which 4.4% comes from the accep- 
tance and operation of the luminosity monitor 
and 4.0% from the uncertainty in the inelastic pp 
cross section. 



II. DETECTOR 

CDF II is a general-purpose detector that is de- 
scribed in detail elsewhere QJ . The components 
relevant to this analysis are briefly described 
here. The detector has an approximate full angu- 
lar coverage with a charged particle tracker inside 
a magnetic solenoid, backed by calorimeters and 
muon detectors. CDF uses a cylindrical coordi- 
nate system in which 9 is the polar angle about 
an axis defined by the proton beam and cj> is the 
azimuthal angle about the beam axis. Particle 
pseudorapidity is defined as r\ = —In tan(0/2). 

The charged-particle tracking system sur- 
rounds the beam pipe and consists of multi- 
ple layers of silicon micro-strip detectors, which 
cover a pseudorapidity region \q\ < 2, and a 3.1 
m long open-cell drift chamber covering the pseu- 
dorapidity region \rj\ < 1 The tracking sys- 
tem is located inside a superconducting solenoid, 
which in turn is surrounded by calorimeters. The 
magnetic field has a strength of 1.4 T and is 
aligned co-axially with the p and p beams. 

The calorimeter system Q is split radially into 
electromagnetic (EM) and hadronic (HAD) sec- 
tions segmented in projective tower geometry, 
and covers the pseudorapidity range \rj\ < 3.6. 
The electromagnetic sampling calorimeters are 
constructed of alternating layers of lead absorber 
and scintillator whereas the hadronic calorime- 
ters use iron absorbers. Proportional cham- 
bers (CES) are embedded in the central elec- 
tromagnetic calorimeter at a depth of about 
6Xo (radiation length), which is the region of 
maximum shower intensity for electrons. In the 
plug region stereo layers of scintillator bars are 
placed at shower maximum. 

A set of drift chambers located outside the cen- 
tral calorimeters (CMU), complemented by an- 
other set behind a 60 cm iron shield (CMP), pro- 
vides muon coverage for \r)\ < 0.6. Additional 
drift chambers and scintillation counters (CMX) 
detect muons in the region 0.6 < |t7|< 1.0 [8j|. 

Multi-cell gas Cerenkov counters [9j located in 
the 3.7 < \r)\ < 4.7 region measure the average 
number of inelastic pp collisions per bunch cross- 



Ill. DATA SAMPLE AND EVENT 
SELECTION 

This analysis is based on an integrated lumi- 
nosity of 2.8 ftr 1 collected with the CDF II de- 
tector between March 2002 and April 2008. The 
data are collected with an inclusive lepton trig- 
ger that requires an electron (muon) with Et > 
18 GeV (pr > 18 GeV/c). The transverse energy 
and transverse momentum are defined as Et = 
E sin 9 and pr = P sin 9, where E is energy mea- 
sured in the calorimeter and p is momentum mea- 
sured by the tracking system. From this inclusive 
lepton data set, events with a reconstructed iso- 
lated electron of Et (muon of pt) greater than 
20 GeV(GeV/c) are selected. Isolation is de- 
fined as the calorimeter energy deposited in a 
cone of radius AR = ^(Arj) 2 + (A</>) 2 = 0.4 
in T) — (f> space centered around the lepton, mi- 
nus the energy deposited by the lepton itself. 
Details on electron and muon identification, or 
lepton ID, criteria used in this analysis are con- 
tained in reference [Til ]. Electrons are identified 
by matching clusters of localized energy deposi- 
tion in the calorimeter to tracks reconstructed 
using hits from the tracking chambers and the 
extra constraint provided by the position of the 
beam line in the transverse direction. We fur- 
ther require that the energy deposition in the 
electromagnetic section of the calorimeter ex- 
ceeds the energy measured in the corresponding 
hadronic section and the lateral cluster energy 
profile agrees with shapes derived from electron 
beam-test data. Muons are identified by match- 
ing tracks to minimum ionizing-likc clusters in 
the calorimeter and to stubs, or sets of radially 
aligned hits, in the muon chambers. Leptons 
passing all of the lepton ID cut and with isolation 
less than 10% of the lepton energy are defined as 
"tight" . They can be of one of four categories: 
electrons reconstructed in the central (CEM) or 
plug (PHX) calorimeter and muons pointing to 
the regions covered by both layers of the two cen- 
tral muon chambers (CMUP) or by the muon 
extension chambers (CMX). These tight leptons 
are also required to be the objects that trigger 
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the event, with the exception of PHX electrons 
which are allowed in events triggered by non iso- 
lated CEM, CMUP or CMX. 

A second electron of Et (muon of pr) greater 
than 20 GeV(GeV/c) is also required using looser 
identification cuts and without the isolation re- 
quirement. "Loose" leptons are either electrons 
or unions which pass the same identification cuts 
as the tight leptons, but fail the isolation require- 
ment. Another category of loose leptons has no 
tight lepton counterpart and is made of muons 
with tracks pointing to regions covered by only 
one of the two central muon chambers (CMU, 
CMP) or with tracks of energy deposition cor- 
responding to a minimum ionizing particle and 
pointing to regions not covered by a muon cham- 
ber (CMIO). CMIO muons, as well as PHX elec- 
trons, must be isolated. 

Each dilepton candidate must contain at least 
one tight lepton and at most one loose lepton. 
These requirements result in 18 different DIL 
dilepton categories, as illustrated in Sec. IIV Al 
where background estimates and acceptances are 
calculated separately for each category. Events 
with more than two tight or loose leptons in the 
final state are rejected as they come mostly from 
background sources like WZ and ZZ events. An- 
other source of trilepton background comes from 
Drell-Yan Z/*y* — > e + e~ events with a radiated 
photon converting into an asymmetric e + e~ pair, 
that is a conversion pair where one electron does 
not reach the minimum track px threshold of 500 
MeV/c needed for the electron to not be trapped 
inside the tracking chamber. The loss in signal 
efficiency from removal of events with three or 
more leptons is only 0.4%. 

A fraction of events passing the dilepton se- 
lection does not originate from pp collisions but 
from beam interactions with the detector and 
shielding material or from cosmic ray sources. 
These events are removed by requiring a recon- 
structed vertex consistent with originating from 
the beam interaction region and within 60 cm of 
the center of the detector along the z-direction. 
We also require that the timing of tracks in 
dimuon events be consistent with both muons 
traveling from the center of the detector outward 
into the tracking chamber [l2| . Electrons from 
conversions of photons in the detector material 
are removed by identifying events with a track 
near the electron track of opposite curvature and 
consistent with coming from a 7 — > e + e~ vertex. 

Jets are reconstructed from the calorimeter 
towers using a cone algorithm with fixed radius 



AR = 0.4 in 77 - <p space [Hj]. The jet E T is 
corrected for detector effects due to calorimeter 
dead zones and to non-linear tower response to 
deposited energy. These effects are convoluted to 
provide the jet energy scale (JES) correction fac- 
tor which estimates the energy of the originating 
parton from the measured energy of a jet 14 1. 



Jets used in the DIL selection are required to 
have corrected Et > 15 GeV and < 2.5. 

We further impose two cuts based on the kine- 
matical properties of the event: the first is a 
cut on the missing transverse energy Ifx 1 which 
measures the transverse energy of the neutrinos 
via the imbalance of the energy detected in the 
calorimeter, after correcting for the presence of 
muons. We require that I$t > 25 GeV, which 
is strengthened to a > 50 GeV requirement if 
any lepton or jet is closer than 20° to the I$t 
direction. This cut, in the following referred to 
as L-cut, is used to reject mainly Z — > t + t~ 
and events with mismeasured Ifr generated by 
jets pointing to cracks in the calorimeter. The 
second, or Z-veto cut, aims at reducing the con- 
tamination of dilepton decays of the Z boson by 
requiring high missing Et significance for ee and 
fifi events with dilepton invariant mass in 76 to 
106 GeV/c 2 region. Missing Et significance, or 
MetSig, is defined as ^ T /\/^| I5r , where Ef? m 
is the sum of transverse energies deposited in 
all calorimeter towers. This variable separates 
events with real I^t due to neutrinos from events 
where the 1$t is due to energy measurement 
fluctations or energy loss in calorimeter cracks. 
These second category of events is expected to 
have a degraded I$t resolution. In the DIL se- 
lection, we use a cut of MetSig > 4 GeV (1/2) . 

Events in the DIL dilepton sample passing the 
L-cut and Z-veto cut become tt candidate events 
if they have at least 2 jets, if the two leptons 
are of opposite charge (OS), and if Ht the trans- 
verse energy sum of leptons, neutrinos and jets is 
greater than 200 GeV. Events in the DIL dilepton 
sample reconstructed with or 1 jet are used as 
control samples for the background estimation. 



1 The scalar quantity Iff is th e magnitude of the missing 
transverse energy vector _9t denned as the opposite of 
the sum over the calorimeter towers J^. E l T rn of the 
transverse energy measured in each tower, Ey, times 
the unit vector in the azimuthal plane that points from 
the beam line to the i th calorimeter tower, n.,. 
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IV. CROSS SECTION 

In this analysis we measure the cross section 
by using the formula 

Nobs — Nbkg /-, N 

atl = AxC (1) 

where N Q bs is the number of dilepton candidate 
events, and Nbkg is the total background. The 
denominator is the product for the acceptance 
for tt candidate events, A, and of the data set 
luminosity L. 

The acceptance, which in our definition in- 
cludes the effects of the detector geometrical ac- 
ceptance, lepton identification and tt to dilepton 
selection, is measured using the pythia Monte 
Carlo program [l5[ to simulate it events of all 
three decay modes (hadronic, lepton + jets, and 
dilepton) with an assumed M top = 175 GeV/c 2 . 
Monte Carlo simulated events are required to 
have both W bosons from top quarks decaying 
to a lepton plus neutrino, where the lepton can 
be either an electron or muon. The tt Monte 
Carlo simulation acceptance prediction of 0.808 
± 0.004 (stat)% is corrected by taking into ac- 
count differences observed between data and the 
Monte Carlo simulation modeling of the detector 
response in independent control samples. The 
following sections describe the implementation 
and checks on the acceptance correction proce- 
dure. 



A. Signal Acceptance 

The available statistics for each sub sample of 
DIL events can be maximized by requiring that 
only detector parts essential for the identification 
of a particular lepton category be fully function- 
ing. For example, PHX electrons require hits in 
the silicon inner vertex detector to reconstruct 
their tracks. Hence the identification of events 
with PHX electrons is limited to data taken with 
"good", i.e. fully functioning, silicon detectors 
but no such requirements is imposed on events 
where the electron is central. To accommodate 
this approach, we rewrite the denominator of the 
cross section formula in equation Q] as the sum 
of the acceptance for each DIL dilepton category 
Ai, weighted by the luminosity relative to that 



category Ci : 

AxC = ^AiXd (2) 

i 

In this analysis four different luminosity sam- 
ples are used, corresponding to the integrated 
luminosity of runs with fully functional sub- 
detectors for trigger CEM electrons and CMUP 
muon, either ignoring the status of the silicon de- 
tectors (2826 pb^ 1 ) or requiring "good" silicon 
(2676 pb _1 ), and runs fully functional also for 
trigger CMX muons, either ignoring (2760 pb _1 ) 
or requiring (2623 pb -1 ) "good" silicon. In defin- 
ing runs good for CEM, CMUP or CMX leptons, 
the distinction between isolated and non-isolated 
leptons is irrelevant. 

The acceptances Ai can be factorized in terms 
of the two leptons l\ and ii comprising the DIL 
category i according to the following formula 

Ai = Ai ± e 2 x Ce 1 e 2 (3) 

where the Ai x i 2 are the raw PYTHIA tt MC ef- 
ficiencies for events with reconstructed leptons 
l\ and £2 passing the full DIL selection and the 
Gi x i 2 are correction factors specific for that lep- 
ton pair. The correction factors are in turn cal- 
culated as: 

= £z x (ftrgi +£trg 2 ~ ^trg 1 ^trg 2 ) X SFf j SFf 2 

(4) 

where e ZQ is an event efficiency while tt rgi and 
SFi i are single lepton trigger efficiency and iden- 
tification efficiency scale factors, respectively. 
The factor e Zo accounts for the efficiency of ± 
60 cm cut on the z-position of the reconstructed 
event vertex. By using a sample of minimally bi- 
ased inelastic interactions, we find that this cut 
accepts 96.63 ± 0.04 (stat)% of the full CDF 
luminous region. The lepton trigger efficiencies 
e tr gi are measured in data samples selected with 
independent sets of triggers and found to be 
around 90% or better. Finally, the scale factors 
SF^ are calculated as ratios of lepton identifica- 
tion efficiencies measured in data and in Monte 
Carlo simulations. 

Table Q] lists all the factors used in the accep- 
tance correction. Their central values are the 
luminosity weighted averages over different data 
taking periods, and the quoted uncertainties are 
only statistical. Using these as inputs to Eq. 
we obtain values ranging from 73% to 93% for the 
correction factors C^ 2 , as shown in Table ILT1 and 
a total denominator ([!]) for the 2.8 fb _1 DIL cross 
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section of 19.43 ± 0.10 pb -1 , where the uncer- 
tainty comes solely from the propagation of the 
statistical uncertainties of each term in Eq. © 
and Eq. ((J). 



B. Check of Acceptance Corrections 

As a cross-check of our acceptance correction 
procedure, we calculate the cross section for Z 
production for each dielectron and dimuon cate- 
gory used in the DIL selection. With this check 
we verify the consistency of the correction proce- 
dure across the different dilepton categories. 

We select events with ee or /i/i in the final 
state and require the two leptons to have oppo- 
site charges and invariant mass in the range 76 
GeV/c 2 to 106 GeV/c 2 . We follow the same lep- 
ton pairing used for the DIL dilepton selection. 
Cosmic ray induced events and events with an 
identified conversion are removed following the 
same criteria used for the tt dilepton selection. 
The number of events selected by these cuts is 
shown in Table [TTTJ 

We use Z/j* — > ee and Z/j* — > pythia 
Monte Carlo simulated samples to calculate the 
raw acceptance of the selection described above 
in the invariant mass 76 GeV/c 2 < Mu < 106 
GeV/c 2 . We use a formulation for the Z cross 
section calculation analogous to the tt cross sec- 
tion formulation in Eq. ([T]). In particular we em- 
ploy the same factorization for the denominator 
correction calculation prescribed in the previous 
section by Eqs. ([2]- [p. 

We perform two checks: time independence of 
the acceptance correction factor for each DIL cat- 
egory and consistency of the correction procedure 
among different categories. With these checks 
we are not trying to measure the Z cross sec- 
tion but rather to determine if our understand- 
ing of the acceptance is correct. The Z cross sec- 
tion has been independently measured by CDF 
as o z = 256 ± 16 pb 0. 

For the first check, we look for possible time 
variations of the measured Z cross section in dif- 
ferent data taking periods corresponding to inte- 
grated luminosities between 200 and 500 pb -1 . 
Figure Q] shows the result of these checks for di- 
electron channels. The error bar in the figure 
reflects uncertainties of data statistics, Monte 
Carlo statistics and the acceptance correction 
procedure. We do not observe any systematic 
trend in the time dependence of the cross section 
for dielectron categories. Same conclusions hold 



for the dimuons channels. As an example, Fig. [2] 
shows where the Z cross section time dependence 
for events with two tight muons. 

For the second check, Table Mil reports the 
cross sections measured over the whole 2.8 fb _1 
data sample for each dilepton category, using the 
efficiencies and scale factors of Table |U The cate- 
gories with two tight leptons (CEM-CEM, CEM- 
PHX, CMUP-CMUP, CMUP-CMX and CMX- 
CMX), which are also the ones with the largest 
acceptance, are consistent with the theoretical 
prediction of 251.6^ pb [l^. Categories with 
a loose muon (CMP, CMU or CMIO) paired to a 
tight muon show some residual variation around 
the average value which is not consistent with 
statistical fluctuations. In order to find a con- 
sistent normalization for all data, we perform a 
fit to the Z cross sections in the different dilep- 
ton categories with three free parameters, cor- 
responding to a multiplicative factor in front of 
the selection efficiency of each of the three loose 
muon categories. The fit returns the cross sec- 
tions reported in the last column of Table Hm and 
an average Z peak cross section of 249.1 ± 0.8 pb, 
as shown in Fig.[3J The Z peak cross section mea- 
sured here does not include the 6% uncertainty in 
the luminosity which is common to all channels 
and which is the dominant uncertainty for the <7z 
measured in [lfj| . The free parameters returned 
by the fit are: SFg MP = 0.977± 0.011, SFg MU = 
1.072±0.011 and SFg MIO = 0.954±0.010. They 
are folded into the acceptance correction proce- 
dure as additional scale factors to be multiplied 
by the lepton identification scale factors of Ta- 
ble U for the appropriate categories. After the 
fit, the single Z cross section measurements are 
consistent with each other within uncertainties, 
with the possible exception of categories contain- 
ing one CMIO loose muon for which we observe 
a maximum deviation equal to 10% of the aver- 
age value. This systematic deviation affects only 
10% of the DIL tt raw acceptance, corresponding 
to the summed contributions of any dilepton pair 
containing a CMIO muon in Table HU Therefore 
we estimate a final 1% systematic uncertainty on 
the acceptance due to the correction procedure. 



V. BACKGROUNDS 

We consider four different sources of standard 
model processes that can mimic the signature 
of dilepton plus 1$t plus 2 or more jets signa- 
ture: diboson events (WW, WZ, ZZ or Wj), 
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TABLE I: Event vertex reconstruction efficiency e zo and list, by lepton type, of trigger efficiency etrg, and 
lepton identification efficiency scale factors SF defined in Eq. Q . These are the luminosity weighted averages 
of efficiencies calculated for different data taking periods. 



Vertex Reconstruction efficiency, e Zo 


0.9663 ±0.0004 


Lepton Type 


Trigger efficiency, e trg 


CEM 


0.965±0.001 


CMUP 


0.917±0.002 


CMX 


0.896±0.002 



Lepton Identification Scale Factor, SF 



CEM 0.987 ± 0.001 

PHX 0.932 ± 0.002 

CMUP 0.927 ± 0.002 

CMX 0.973 ± 0.002 
CMU 0.97 ± 0.01 

CMP 0.90 ± 0.01 

CMIO 0.99 ± 0.01 



TABLE II: List, by dilepton category, of raw acceptance Ae 1 e 2 , correction factor Ci t i 2 and luminosity d 
used in the calculation of the denominator for the 2.8 fb _1 DIL cross section measurement. The acceptance of 
each category includes contributions from non isolated loose leptons. The Ae 1 e 2 uncertainty comes only from 
the MC statistics. The error in the Ct 1 t 2 comes from the propagation of the dilepton efficiency uncertainties 
of Table U 



DIL Category A tl t 2 {%) C tl e 2 gggb^J 



CEM-CEM 





1224 


± 





.0017 





.9338 


zb 





.0019 


2826 


CEM-PHX 





.0470 


± 





.0010 





.8658 


dz 





.0027 


2676 


CMUP-CMUP 0. 


.0498 


± 





.0011 





.8189 


dz 





.0025 


2826 


CMUP-CMU 





.0191 


± 





.0007 





.7920 


dz 





.0040 


2826 


CMUP-CMP 





.0267 


± 





.0008 





.7299 


dz 





.0035 


2826 


CMUP-CMX 





.0474 


± 





.0010 





.8569 


zb 





.0020 


2760 


CMUP-CMIO 





.0234 


± 





.0007 





.8050 


zb 





.0040 


2826 


CMX-CMX 





.0106 


± 





.0005 





.8996 


zb 





.0027 


2760 


CMX-CMU 





.0075 


zb 





.0004 





.8134 


zb 





.0043 


2760 


CMX-CMP 





.0115 


zb 





.0005 





.7495 


zb 





.0037 


2760 


CMX-CMIO 





.0101 


± 





.0005 





.8267 


zb 





.0043 


2760 


CEM-CMUP 





.1769 


± 





.0020 





.8737 


zb 





.0018 


2826 


CEM-CMU 





.0349 


± 





.0009 





.8879 


zb 





.0040 


2826 


CEM-CMP 





.0475 


± 





.0010 





.8182 


zb 





.0035 


2826 


CEM-CMX 





.0845 


dz 





.0014 





.9171 


zb 





.0019 


2760 


CEM-CMIO 





.0410 


dz 





.0010 





.9025 


zb 





.0041 


2826 


PHX-CMUP 





.0327 


dz 





.0009 





.7723 


zb 





.0029 


2676 


PHX-CMX 





.0147 


± 





.0006 





.7922 


zb 





.0032 


2623 
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Dataset Periods 



FIG. 1: Z cross section using dilcpton categories with two electrons as a function of data taking period. 




Dataset Periods 



FIG. 2: Z cross section using dilepton categories with tight muons as a function of data taking period. 
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FIG. 3: Z peak cross section measured for each inclusive same flavor dilepton category using the full 2.8 fb~ 
data sample after applying the loose muon scale factor. Not included is a 6% uncertainty in the luminosity 
measurement common to all channels. The band represents the mean value of 249.1 ± 0.8 pb. 



TABLE III: Number of selected events and Z cross section in 2.8 fb -1 for the different ee and fifi dilepton 
categories with dilepton invariant mass in the range 76 GeV/c 2 to 106 GeV/c 2 . Results are given both for the 
original and the fitted Z cross sections. The cross section uncertainties are only from the data statistics and 
from the propagation of the uncertainty in the single lepton efficiency of Table [I] 



Category 


Number of Events Z Cross Section (pb) Orh; 


>inal Fitted 


CEM-CEM 


50519 


248.0 ± 1.4 


248.0 ± 1.4 


CEM-PHX 


51468 


250.3 ± 1.4 


250.3 ± 1.4 


CMUP-CMUP 


14096 


251.2 ± 2.4 


251.2 ± 2.4 


CMUP-CMU 


5914 


261.4 ± 3.7 


243.8 ± 4.5 


CMUP-CMP 


6944 


244.3 ± 3.2 


250.1 ± 4.2 


CMUP-CMX 


20007 


247.1 ± 2.0 


247.1 ± 2.0 


CMUP-CMIO 


9099 


229.2 ± 2.7 


240.3 ± 3.6 


CMX-CMX 


5563 


250.8 ± 3.6 


250.8 ± 3.6 


CMX-CMU 


3969 


275.5 ± 4.7 


257.0 ± 5.3 


CMX-CMP 


4280 


241.8 ± 3.9 


247.6 ± 4.8 


CMX-CMIO 


5771 


251.0 ± 3.6 


263.2 ± 4.4 



Drcll-Yan production of tau leptons (DY— » tt), 
Drcll-Yan production of electrons or muons with 
additional Ifx (if the event is an actual Drell- 
Yan event, there is no 1$t so we refer to this as 
fake Ifx ) (DY— > eej \x\x) and QCD production of 
W boson with multiple jets in which one jet is 
misidentified as a lepton ('W+jet fakes"). The 
two dominant sources of background are DY— > 
ee/n/i and VF+jet fakes. These two processes 
have production cross sections much larger than 
the ti, but they can only contaminate the tt dilep- 
ton signature of two leptons plus jets and large 



Ijj r p when misreconstructions of the event create 
either some large fake or a jet misidentified as 
a lepton. Because it is difficult to use the Monte 
Carlo simulation to predict the effect of event 
misreconstruction in our detector, we estimate 
the background from these two processes using 
data-based methods, as discussed in Sec. IV Al and 
IV Bl respectively. The diboson and DY— > tt 
backgrounds are calculated using Monte Carlo 
simulation expectation as described in Sees. IV CI 
and IV D1 Corrections are applied for trigger and 
lepton ID efficiencies following the same proce- 
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dure described in Sec. II V A I 

Our strategy for validating the background es- 
timation is to compare data and background es- 
timates in the 0-jet and 1-jet bins, as discussed 
in Sec. IVlTl 



A. VK+jet fakes 

Events with a single W boson plus jets can 
simulate the dilepton signature when one of the 
jets is misidentified as a lepton. The W+jet fake 
contamination is calculated in two steps: first we 
extract the probability of a generic QCD jets fak- 
ing the signatures of different lepton categories; 
then we apply these probabilities to weight events 
in the data containing one and only one high pt 
lepton plus jets. 

The fake probabilities are measured in generic 
jets from QCD decays by selecting "fakeable" lep- 
tons, which are jets passing minimal lepton iden- 
tification criteria described below. We do not 
consider separately the heavy flavor contribution 
to our backgrounds because the probability for a 
b or c quark to become a well reconstructed high 
Pt lepton is very small. We define different cate- 
gories of fakeable lepton, one per high px lepton 
category in the DIL dilepton selection. 

Jets with a large fraction of neutral to charged 
pion production can create signatures with low 
track multiplicity and large energy deposition in 
the electromagnetic calorimeter, thus faking the 
presence of electrons. We define fakeable elec- 
trons as tracks of pr > 20 GeV/c pointing to an 
electron-like cluster with energy deposition in the 
electromagnetic section of the calorimeter far ex- 
ceeding the energy measured in the hadronic sec- 
tion, namely with Ehad/Eem < 0.125. Fakeable 
electrons are further divided into objects that 
can fake CEM or PHX electrons depending on 
whether their clusters belong to the central or 
plug section of the calorimeter. We label them 
TCEM and TPHX, respectively. Fakeable for the 
non isolated electrons do not require isolation for 
the central cluster and are called NCEM. 

Jets whose full hadronic activity is limited to 
single charged pions or kaons with a late shower 
development or decay in flight might deposit lit- 
tle energy in the calorimeter but generate hits 
in the muon chambers, thus faking the signature 
of a muon. We define fakeable muons as good 
quality tracks of pr > 20 GeV/c with E/p < 1. 
Depending on which muon sub-detectors these 
tracks point to, we label as TCMUP, TCMX, 



LMIO and LMUO fakeable muons which can 
fake tight CMUP, tight CMX, loose CMIO or 
loose CMU/CMP muons, respectively. Fakeable 
muons that fail the isolation requirement are put 
together into a single NMUO category as long as 
they point to any muon sub-detector. 

We select fakeable leptons among generic QCD 
jets collected in four different control samples, 
whose main trigger requirement is the presence 
of at least one jet of E^ rg > 20, 50, 70 and 100 
GeV, respectively. The simplification of the jet 
algorithm used in these trigger selection tends 
to underestimate the energy of the offline recon- 
structed jets. To ensure a trigger efficiency of 
90% or greater we require the trigger jet to have 
reconstructed Et greater than 35, 55, 75 and 105 
GeV, respectively, in the four jet samples. The 
resulting probabilities are labeled Jet20, Jet50, 
Jet70 and JetlOO fake lepton probabilities. To 
minimize real lepton contamination, we require 
that fakeables in the denominator of the fake 
probability fail one or more of the standard lep- 
ton identification cuts. For the numerator in- 
stead we require that the fakeable leptons pass 
all of the lepton identification requirements. We 
estimate the contamination of real leptons from 
Ws, Drell-Yan or dibosons, using Monte Carlo 
simulation predictions for the number of events 
with one lepton and at least one jet above the 
Et threshold. 

We use the fake lepton probability measured 
in the Jet50 sample as our primary estimator 
to apply to data events because the jet energy 
spectrum in the Jet50 sample is the closest to 
the energy spectrum of jets in the dilepton plus 
missing Et sample. The fake probabilities for 
different lepton categories show a dependence on 
the transverse energy of the fakeable lepton. To 
properly account for difference in the pt spec- 
trum of fakeable leptons in QCD jets vs W+jets, 
we calculate fake probabilities in six px ranges as 
shown in Table ITV1 

The uncertainties on the fake probabilities in 
Table IIVI are only statistical. Variations in fake 
probabilities between the different QCD jet sam- 
ples are used to estimate a systematic uncertainty 
in the lepton fake estimate. Figure 0] shows a 
comparison between the number of fake lepton 
events observed in the Jet20, Jet70 and JetlOO 
data sample, after integration over the full j>t 
spectrum, and the number predicted by the Jct50 
fake probabilities of Table IIVI We assess a 30% 
systematic uncertainty on the ability of the Jet50 
fake probabilities to predict electron and muon 
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fake contamination in samples with a wide range 
of jet energy. 

We define "lepton+fakeable" as those events 
in the central high px lepton data sets with one 
and only one good high px lepton, tfx >25 GeV 
and a second fakcable object failing at least one 
standard lepton identification cuts. The fakeable 
object, which can be from any of the fakeable 
categories defined above, is paired to the good 
lepton and treated as the second lepton in the 
event when calculating any of the kinematic vari- 
ables used in the top quark DIL selection, such as 
dilepton invariant mass, corrected I^t, and Ht- 
Jets found in a cone of AR < 0.4 around the fake- 
able lepton are not included in the jet multiplicity 
count of that event because those jets are associ- 
ated with the fake lepton in this VF+jet fake esti- 
mation scheme. The fake lepton contamination is 
calculated by weighting each "lepton+fakeable" 
event found in data by the fake probability in 
Table IIVI If more than one fakeable object is 
found in the event, we pair each of them to the 
good lepton and add their single fake contribu- 
tions. The fake dilepton background thus calcu- 
lated contains a statistical component, which is 
the sum of the fake probability uncertainty itself 
and the statistics of the "lepton+fakeable" sam- 
ple. 

As a check, we compare the same sign 
"lepton+fakeable" prediction to the number of 
VF+jet fakes with same sign dilepton candidates 
in the signal regions. We define as fake lepton 
charge the charge of the track associated to the 
fakeable lepton. Same sign dilepton candidates 
are corrected for the presence of same sign pairs 
coming from it, DY or diboson events that are 
simulated in our Monte Carlo simulations. The 
results of this check are shown in Table [Vj Al- 
though the fi/i channel shows deviation at the 3 
standard deviation level for some jet multiplicity 
bins, the agreement in the final predictions over 
all dilepton categories is at the one standard de- 
viation level. 



B. Drell-Yan to ee/nn Background 

The contamination from Z/j* ee/fifi de- 
cays is calculated using a combination of data 
based and MC based predictions. We define DIL 
data samples enriched in DY events after the L- 
cut by inverting the Z-veto cut and extrapolat- 
ing the remaining DY contamination in the sig- 
nal region by using the relative contribution of 



Z/j* — > ee/ nn decays passing and failing the Z- 
veto cut as predicted from MC. The Z-veto cut 
(see Sec. IIIII) requires that the dilepton invariant 
mass be outside the Z window region of 76 to 106 
GeV/c 2 , or, if inside, that the event have missing 
E T significance tyr/y^E^ > 4 GeV (1 / 2 \ 

We calculate the DY— >ee//z/u contamination as 
the sum of two contributions, one outside the Z 
window region, N out , and one inside the Z win- 
dow with high MctSig, Nhigh- The first contribu- 
tion is calculated as: 

Nout = Rout/i„ (N° T - Nf n KG ) (5) 

where N? T and N? KG represent the number of 
events inside the Z-window passing the L-cut 
in data and in non-DY MC background predic- 
tions, respectively. R ou t/m is the ratio of Z /^* — > 
ee / /Lt/Lt events outside to inside the Z window pre- 
dicted by the ALPGEN [lj| Monte Carlo genera- 
tor. 

The second contributions is calculated as: 

Nhigh = Rhigh/low( N !ow _ N k,w ) (6) 

where and N^ G represent the events in- 

side the ^-window passing the L-cut with Met- 
Sig < 4 GeV^ 1 / 2 ) for data and for non-DY MC 
background predictions, respectively. Rhi g h/iow is 
the ratio of events passing/failing the MctSig > 
4 GeV (1 / 2 ) cut predicted by ALPGEN. Table [VTl 
summarizes the inputs to Eqs. ([5]) and ([6]) and the 
final values of N ou t and Nhigh for each jet multi- 
plicity bin. For the calculation of it contribution 
to N BKG we use the prediction of 6.7 pb for the 
cross section. We later correct this iteratively to 
the value measured in the data. 

The DY contamination in the signal sample is 
extracted from the N out and Nhigh estimates in 
the > 2 jet bin, corrected for the efficiency of the 
Ht > 200 GeV and of the opposite sign lepton 
cuts. The combined efficiency for these two cuts 
is calculated using ALPGEN simulated Z samples 
and shown as eos in Table I VIII 

The contamination of Z/j* — > /i/i to e/j, events 
comes mostly from cases where one of the fi- 
nal state muon radiates a very energetic photon. 
These photons, which are almost collinear to the 
muon, deposit their energy in the EM calorime- 
ter and produce a cluster which is associated to 
the original muon track and fakes the electron 
signature. The missed muon gives rise to a siz- 
able Iff in the event, curtailing the effectiveness 
of the L-cut and Z-veto to reject them. As no 
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TABLE IV: Jet50 fake probabilities vs fakeable lepton pr for different fakeable categories. The uncertainties 
are statistical only. Due to the definition of fake probability, the denominator can fluctuate to be smaller than 
the numerator in low statistics high pr bins, hence fake rate values exceeding 100%. 



Jet50 Fake Probabilities (%) in px range (GeV/c) 

Fakeable [20-30] [30-40] [40-60] [60-100] [100-200] > 200 

TCEM 4.97 ± 0.09 3.68 ± 0.08 2.39 ± 0.01 2.88 ± 0.02 3.26 ± 0.11 4.98 ± 3.44 

NCEM 0.74 ± 0.1 0.53 ± 0.01 0.58 ± 0.01 0.49 ± 0.01 0.80 ± 0.07 

TPHX 12.6 ± 0.1 13.9 ± 0.1 12.7 ± 0.1 20.4 ± 0.2 26.8 ± 0.2 65.0 ± 57.3 

TCMUP 0.99 ± 0.04 2.30 ± 0.10 3.82 ± 0.14 6.30 ± 0.28 4.13 ± 0.47 

TCMX 0.91 ± 0.06 2.20 ± 0.15 4.74 ± 0.21 5.09 ± 0.55 0.76 ± 1.98 0.44 ± 0.25 

LMUO 2.48 ± 0.05 2.88 ± 0.14 4.41 ± 0.23 5.45 ± 0.47 7.51 ± 0.64 0.41 ± 0.26 

LMIO 21.0 ± 0.1 25.4 ± 0.3 27.8 ± 0.4 40.1 ± 0.1 37.0 ± 1.5 107.1 ± 14.8 

NMUO 0.51 ± 0.01 0.36 ± 0.01 0.23 ± 0.01 0.25 ± 0.01 0.29 ± 0.04 0.17 ± 0.18 



TABLE V: Comparison between the same sign dilepton fake background prediction using the fake rate tables 
and the numbers of same sign dilepton candidates found in the signal region, after MC subtraction of standard 
model contamination sources. 



Number of SS Dilepton Events ee fj,fi e/j, all 



Corrected Candidates 


in 





jet 


13.0 


± 


4.4 


11.3 


± 


3.7 


14.8 


± 


5.0 


39.0 


± 


7.6 


Predicted Candidates 


in 





jet 


8.1 


± 


2.6 


7.8 


± 


2.7 


18.0 


± 


5.6 


33.9 


± 


10.4 


Corrected Candidates 


in 


1 


jet 


7.8 


± 


3.5 


0.5 


± 


1.4 


24.9 


± 


5.6 


33.2 


± 


6.7 


Predicted Candidates 


in 


1 


jet 


5.1 


± 


1.6 


6.9 


± 


2.2 


25.4 


± 


7.8 


37.3 


± 


11.4 


Corrected Candidates 


in 


2 


or more jets 


5.0 


± 


2.7 


0.7 


± 


1.0 


24.7 


± 


5.2 


30.4 


± 


5.9 


Predicted Candidates 


in 


2 


or more jets 


3.8 


± 


1.2 


7.5 


± 


2.4 


24.5 


± 


7.6 


35.9 


± 


11.0 



data based control sample is available for this 
contamination, we estimate it using Monte Carlo 
simulation predictions. 



C. Diboson Background 

The diboson processes, WW, WZ, ZZ and 
Wj, can mimic the signature of the tt signal via 
different mechanisms, with real leptons and I^t 
from W and Z decays and jets produced by boson 
hadronic decays or initial and final state radia- 
tion. For WW events, the two leptons and the 
l^T are produced when both W's decay semi- 
lcptonically but the jets require some hadronic 
radiation external to the diboson system. For 
WZ and ZZ events, the two leptons come from 
the Z boson while the other W or Z boson pro- 
vides the jets via their hadronic decays. As these 
decays do not contain any neutrino, some mech- 
anism to produce fake Jf*T is required. Finally 
for W^f events, one lepton plus Ifx is generated 
from the semi-leptonc W decay while the second 
lepton is produced from an asymmetric 7 conver- 
sion in which one of the two electrons has little 



energy and is caught spiralling inside the central 
drift chamber. Like in the WW case, the Wj 
system is accompanied by hadronic jets. Events 
involving VF+jets fake leptons, with a real lepton 
from W boson paired to a fake lepton from the 
hadronic decays of the other boson, are removed 
from the MC to avoid double counting. 

Only WW background contribute to the ee, 
/i/j, and e/i final states in the same proportion as 
the ti signal. Diboson processes involving a Z 
contribute preferentially to the same flavor lep- 
ton categories. W-f events do not contribute any 
background to the /i + [iT category given the neg- 
ligible probability that the photon will convert to 
a muon pair. 

The WW, WZ and ZZ processes are simu- 
lated with the pythia Monte Carlo generator. 
Their production cross section is taken from the 
latest next-to-leading order (NLO) MCFM ver- 
sion [17] and CTEQ6 dS] PDF predictions to be 
cr ww = 12 A ± 0.8 pb, a wz = 3.7 ± 0.1 pb. For 
the ZZ events, a cross section ozz = 3.8 pb is 
assumed with an uncertainty of 20%. W"/ decays 
are simulated with the BAUR Monte Carlo gen- 
erator [19(. The leading order (LO) production 
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Using jet50 fake rates 
MET>25 GeV lepton + fakeable events 



TCEM 
NCEM 
TPHX 
TCMUP 
TCMX 
LMUO 
LMIO 
NMUO 
jet20: total 

TCEM 
NCEM 
TPHX 
TCMUP 
TCMX 
LMUO 
LMIO 
NMUO 
jet70: total 

TCEM 
NCEM 
TPHX 
TCMUP 
TCMX 
LMUO 
LMIO 
NMUO 
jetl 00: total 



-0.079 +0.036 

0.213+0.007 

-0.147+0.014 

-0.102+0.134 

-0.374+0.308 

-0.663 ±0.377 

0.415+0.013 

-0.181+0.015 

0.076+0.022 

0.375+0.008 
-0.091+0.001 
0.198+0.003 
-0.017+0.030 
0.104+0.111 
-0.069 ±0.037 
-0.010 ±0.008 
0.1 37 ±0.001 
0.118±0.011 

0.626 ±0.004 
-0.126 ±0.001 
0.441 ±0.002 
-0.109 ±0.032 
0.049 ±0.1 16 
0.047 ±0.031 
-0.1 82 ±0.009 
0.389 ±0.001 
0.273 ±0.008 



-1 1 

(Obs-Pred)/Obs 



FIG. 4: Ratio of observed total number of fake leptons for each fakeable category vs the Jet50-based prediction 
normalized by the number of observed. The predictability of the jet50 Pt dependent fake rate is good at the 
30% level, as shown by the band in the plot. When error bars are not shown they are smaller than the dot 
size. 



cross section of crw-y = 32 ± 3 pb is assumed, and 
multiplied by a K-factor of 1.36 20] to correct for 
NLO effects. The Wj Monte Carlo generator ac- 
ceptance prediction is multiplied by a conversion 



inefficiency scale factor of 1.15 ± 0.35 to correct 
for the imperfect simulation of the tracking vari- 
ables used in the conversion identification algo- 
rithm. 



16 



TABLE VI: Inputs to Eqs. ([5} and (0 and for each dilepton flavor and jet multiplicity. N out and N n i g h are 
the final values of the DY— > ee and fj,fi background contamination outside the Z peak region and inside the Z 
peak region with high MetSig, respectively. 





0-jet 
ee im 


1-jet 
ee /i/i 


>2-jets 
ee fifi 


N BKG 
in 

R'out /in 


78 

29.6±1.4 
0.39±0.05 


45 
20.1±1.0 
0.45±0.09 


76 
16.5±1.6 
0.31±0.05 


58 

11.9±1.4 
0.31±0.06 


73 

17.0±1.2 
0.32±0.02 


43 
16.8±1.0 
0.26±0.02 


Nout 


19.1±3.7 


11.3±3.2 


18.2±2.8 


14.2±2.5 


18.2±2.8 


6.8±1.7 


low 
-^high/low 


65 
12.6±0.7 
0.026±0.009 


37 
9.5±0.6 
0.010±0.005 


69 
5.9±0.5 
0.049±0.006 


54 
5.5±0.5 
0.049±0.011 


65 
6.7±0.2 
0.040±0.003 


39 
8.0±0.5 
0.040±0.004 


Nhigh 


1.23±0.21 


0.26±0.06 


2.85±0.41 


2.34±0.36 


2.16±0.32 


1.22±0.25 



TABLE VII: Ht and opposite sign cut efficiency for the DY— > ee and fifi background contamination in > 2 
jet region. The efficiency is calculated separately for events outside the Z peak region passing the L-cut, and 
for events inside the Z peak region also passing the MetSig > 4 GeV (1/2) cut. 



£h t ,os ee MM 

for N out events 0.54±0.03 0.60±0.05 
for N in events 0.95±0.01 0.99±0.01 



Monte Carlo generators do not correctly model 
the jet production from hadronic radiation, as is 
seen by comparing the jet multiplicity spectra of 
data and MC predictions for ee and events in 
the Z peak region. Data, even after correcting 
the jet multiplicity spectrum for other SM con- 
tributions, have higher fractions of events in the 
2 or more jet bins compared to predictions. We 
calculate jet multiplicity scale factors CNj as ra- 
tios of data and MC events in each jet bin, after 
normalizing the MC to the number of data in the 
Z peak region. These scale factors, shown in Ta- 
ble IVIIIl are used to correct the jet multiplicity 
of WW and Wj events. A 5% systematic uncer- 
tainty on this correction is assessed by comparing 
jet multiplicity scale factors calculated with dif- 
ferent generators. 

D. Drell-Yan — > tt Background 

Z/j* t + t~ decays are simulated with the 
ALPGEN generator. These events can fake the 
dilepton plus Jf,T plus 2 or more jets signa- 
ture when both t's decay semi-leptonically to 
£ + V£i7 T £~ viv T and jets from initial and final state 
radiation are present. The contamination from 



TABLE VIII: Jet multiplicity scale factors for Z — s> 
e + e~ and Z — \ pT events in the 0-jet bin (Coj), 
1-jet bin (Cij), and > 2-jet bin (C2j), respectively. 
The last column is the weighted average of the two 
same flavor Z samples and it is used as the correction 
factor for e/x reconstructed events. The uncertainties 
shown here are statistical only. 



Jet multiplicity Scale Factor 
e+e" /x+M" 
C oj 1.017 ± 0.010 0.999 ± 0.011 1.010 ± 0.010 
Cij 0.918 ± 0.012 0.991 ± 0.012 0.948 ± 0.008 
C 2i 1.056 ± 0.020 1.123 ± 0.020 1.082 ± 0.014 



this process is expected to contribute equally to 
the e + e~ and /J. + fi~ categories and to be twice 
as big in the e^ channel. The neutrinos from 
the semi-leptonic r decays tend to have lower en- 
ergy than the neutrinos in the tt dilepton sample 
and align along the direction of the leptonic de- 
cay when the Z recoils against the external jets. 
Hence a big fraction of the Z/'y* — > tt events 
are removed by the L-cut, the cut on the event 
Ifx > 25 GeV or I$t > 50 GeV in case any lepton 
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or jet is closer than 20° to the Ifr direction (see 

secnm. 

The final contamination from this process is es- 
timated using a Monte Carlo simulation and as- 
sumes a Z — y tt cross section of 251. 6jl§'f pb |10j . 
The simulated samples are generated using ALP- 
GEN generator [l|| that has built-in matching of 
the number of jets, coupled with PYTHIA |l5j for 
the shower evolution and evtgen [21| for the 
heavy- flavor hadron decays. All simulated events 
were run through the full CDF detector simula- 
tion. To correct for NLO effects, this value is 
further multiplied by a K- factor of 1.4 [22j . The 
MC predictions in the different jet bins arc finally 
rescaled by the Cisrj scale factors, as discussed in 
Sec. V C. 



VI. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainty for the cross sec- 
tion measurement has two main contributions: 
systematics in the tt dilepton acceptance and sys- 
tematics in the background estimation. We dis- 
tinguish between the uncertainties affecting only 
the signal or the background from the uncertain- 
ties common to both. 

For the signal acceptance, we consider sys- 
tematic uncertainties coming from different MC 
generators, different assumed amounts of ini- 
tial (ISR) and final (FSR) state radiation in 
the PYTHIA Monte Carlo, calculated by compar- 
ing data to expectations for the pr spectrum 
of the dilepton system in Drell-Yan events and 
for the kinematic distributions of the underlying 
events and different parton distribution functions 
(PDF). These sources are uncorrelatcd from each 
other. The two remaining and largest sources of 
acceptance systematics are common to signal and 
background Monte Carlo simulation predictions. 
They arise from uncertainties in the lepton iden- 
tification (ID) scale factors and jet energy scale 
(JES). Comparing the lepton ID scale factors cal- 
culated for Z events with 0, 1 and > 2 jets, we de- 
rive a systematic associated to the Monte Carlo 
generator acceptance correction of 2%. This is 
added in quadrature to the 1% systematic un- 
certainty on the acceptance correction procedure 
derived from measurement of the Z cross section 
in different dilepton channels (see Sec. II V A"]) , for 
a total systematics on the lepton ID correction of 
2.2%. The JES uncertainty is calculated by mea- 
suring the shift in acceptance due to varying the 
jet energy scale correction applied to each jet in 



the event by ±1 standard deviation of its system- 
atic uncertainty. Table IIXI summarizes the sys- 
tematic uncertainties in the tt acceptance sepa- 
rated by the contributions which are independent 
and contributions which are common with the 
systematic uncertainty in the background pre- 
diction. Common contributions affect both the 
numerator and the denominator of Eq. ([1]) used 
to calculate the final tt cross section. Their corre- 
lation is taken into account when calculating the 
systematic uncertainty on the measured tt cross 
section. 

TABLE IX: Summary of systematic uncertainties af- 
fecting the tt acceptance: the entries in the top part 
of the table are treated as uncorrelated and added in 
quadrature when calculating their contribution to tt 
cross section systematic uncertainty via the denomi- 
nator of Eq. {1} ; the two entries in the bottom part of 
the table are also added in quadrature but their cor- 
relation with the systematic uncertainty in the back- 
ground prediction, which appears in the numerator 
of Eq. |T]|, is taken into account when calculating the 
final systematic uncertainty in the ti cross section. 



Source 


Systematic Error (%) 


MC Generator 


1.5 


ISR 


1.7 


FSR 


1.1 


PDF 


0.8 


Lepton ID 


2.2 


Jet Corrections 


3.2 



Uncorrelated sources of systematic uncertain- 
ties affecting the backgrounds are the 30% sys- 
tematic uncertainty on the fakes contamination, 
the 30% uncertainty on the conversion ineffi- 
ciency scale factor affecting the W"f contami- 
nation and the theoretical uncertainties, rang- 
ing from 2% to 10%, on the production cross 
sections of diboson and Z— > tt processes. Al- 
though large, each of these systematics affects 
only a fraction of the total background. Finally, 
a systematics common to most Monte Carlo gen- 
erator predictions of background processes with 
jet production from QCD radiation comes from 
the 5% uncertainty in the Cnj correction factors 
of Sec. EH] 

VII. CONTROL SAMPLES 

We use dilepton events passing both the L-cut 
and Z-veto cut, but with only or 1 jets as con- 
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trol samples for the background calculation. The 
ti contamination to the jet samples is negligible 
and the contribution to the 1 jet sample is small. 
The results for the number of observed and ex- 
pected events in the various dilepton flavor cate- 
gories are shown in Table [X] The ti contribution 
has been calculated assuming a production cross 
section of 6.7 pb. For the 0-jet bin data, the 
single largest contribution is from diboson pro- 
duction, followed by W -1-jet fakes and DY pro- 
duction; for the 1-jet bin data the single largest 
contribution is from W+jet fakes, followed by di- 
boson and DY production. Figures [5] and H] show 
the 0-jet and 1-jet data overlaid on top of the 
background and tt predictions for four different 
kinematic distributions: single lcpton px, dilep- 
ton invariant mass, event missing transverse en- 
ergy Iff and total scalar transverse energy H/p. 
The largest deviation, still below the 2 standard 
deviation level, is in the Nj et = control sample 
for the [i/i channel. Overall the data are in good 
agreement with the expectations of background 
plus signal (tt). The agreement is quantified in 
terms of the probability for the x 2 1 ndf distribu- 
tion and shown as "x 2 Test" on the figures. 

VIII. tt USING 2 JET SELECTION 

As an intermediate step toward the final re- 
sult, Table IXII shows the predictions for signal 
and background in events with two or more jets 
passing all of the tt selection criteria, but the 
Ht > 200 GeV and the opposite lepton charge 
requirement. For this sample the tt signal con- 
tribution is almost equal to the total background 
contribution. There is good agreement between 
predicted and observed number of events both in 
overall normalization and in the bin-by-bin dis- 
tribution for the same four kinematic variables 
used in the and 1-jet control samples, as shown 
in Fig. [7] The agreement is quantified in terms 
of the probability for the x 2 /ndf distribution and 
shown as "x 2 Test" on the figures. 



IX. RESULTS 

The signal and background DIL candidate 
events, that is events in the 2 or more jet sam- 
ples after the final Ht > 200 GeV and oppo- 
site charge requirements, are shown in Table IXTfl 
separately for the different dilepton flavor contri- 
bution. The tt rate is computed assuming a ti 



production cross section in agreement with the 
NLO standard model calculation for a top mass 
of 175 GeV/c 2 , of 6.7±££ pb 0. The sum of the 
background and signal contributions is labelled 
"Total SM expectation" and can be compared 
to the number of "Observed" data in 2.8 fb _1 . 
Fig. [5] shows the ti and the different backgrounds 
overlaid on the data, for the single lepton pt, 
dilepton invariant mass, the event I$t and Ht 
distributions. Again there is overall good agree- 
ment between data and total background plus ti 
expectations, as shown by the probability for the 
X 2 /ndf distribution reported on the figures. 

Table IXIIII summarizes the background and 
signal predictions for the 0, 1 and > 2 jet-bin 
control samples and for the signal sample. Figure 
IH1 shows the overall number of candidate events 
in the different jet multiplicity bins overlaid on 
top of a stacked histogram of the different back- 
ground components. The band gives the ti con- 
tribution for a cross section of 6.7 pb. The 
hatched area represents the uncertainty in the 
total background estimate. From the difference 
between the observed data and the total back- 
ground predictions, we measure a cross section 
for ti events in the dilepton channel of: 

ait = 6-27 ± 0.73(stat) ± 0.63(syst) ± 
0.39(lum) pb. 

where the first uncertainty is statistical, the sec- 
ond is the convolution of the acceptance and 
background systematics and the third comes 
from the 6% uncertainty in the luminosity mea- 
surement. This results assumes a top quark mass 
of M t = 175 GeV/c 2 . Studies of the DIL se- 
lection acceptance vs Mj shows an increase in 
acceptance of 3% for each 1 GeV/c 2 , in the 
range ± 2 GeV/c 2 around the combined Tevatron 
top quark mass measurement of M t = 173.1 ± 
0.6 s tat±l-l sys t GeV/c 2 [23[. The theory cross sec- 
tion decreases by approximately 0.2 pb for each 
1 GeV/c 2 increase over the mass range from 170 
to 180 GeV/c 2 . 

As a test of lepton universality, we quote the 
results for the individual dilepton flavor decay 
modes: 

(T t t{ee) = 4.57 ± 1.56(stat) ± 0.58(syst) pb 
0ti(W*) = 7 - 47 ± l-63(stat) ± 0.79(syst) pb 
a t t(en) = 6.43 ± 0.95(stat) ± 0.69(syst) pb. 

All of the results are consistent with each other. 
Similar conclusions hold for the cross section of 
signal events where both leptons are isolated, 
which is a sample extensively used in other SM 
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TABLE X: Summary table, by lepton flavor content, of background estimates, tt predictions and observed 
events in data corresponding to an integrated luminosity of 2.8 fb for the jet (top) and 1 jet (bottom) 
bins, respectively. The quoted uncertainties are the sum of the statistical and systematic uncertainty. 



Nj et = Control Sample per Dilepton Flavor Category 



Source 


ee 


fifj, 


e/j, 


It 


WW 


36.36±3.26 


30.10±2.71 


76.40±6.76 


142.87±12.57 


wz 


2.88±0.21 


4.56±0.31 


4.21±0.29 


11.65±0.74 


zz 


4.13±3.19 


4.25±3.28 


0.41±0.32 


8.79±6.78 




14.26±4.90 


O.OOiO.OO 


13.77±2.63 


28.03±7.11 


DY^ TT 


0.95±0.26 


0.79±0.23 


2.15±0.39 


3.89±0.58 


DY^ ee + fi/i 


20.32±3.94 


11.59±3.34 


8.13±1.33 


40.04±5.77 


W+jet fakes 


18.72±5.73 


15.29±4.93 


38.47±11.72 


72.49±19.02 


Total background 97.64±14.47 66.58±9.23 143.54±15.65 307.76±35.84 


tt (a = 6.7 pb) 


0.15±0.03 


0.18±0.03 


0.34±0.04 


0.67±0.06 


Observed 


99 


96 


147 


342 



Njet = 1 Control Sample per Dilepton Flavor Category 



Source 


ee 


nn 


e/i 


It 


WW 


9.74±1.08 


8.73±0.97 


21.55±2.31 


40.01±4.24 


WZ 


4.95±0.25 


2.66±0.15 


4.11±0.21 


11.72±0.52 


ZZ 


1.59±1.23 


1.58±1.22 


0.91±0.70 


4.08±3.14 


Wj 


3.70±1.47 


O.OOiO.OO 


4.43±1.12 


8.14±2.27 


DY^ TT 


4.64±0.83 


4.42±0.78 


8.81±1.50 


17.87±2.99 


DY^ ee + l-ifJ, 


21.04±4.27 


16.48±3.56 


3.31±0.90 


40.83±7.28 


VF+jct fakes 


12.14±3.73 


14.84±4.60 


67.26±20.43 


94.23±26.16 


Total background 


57.80±8.97 


48.70±7.57 110.37±21.27 216.87±32.46 


tt (a = 6.7 pb) 


3.94±0.22 


4.02±0.22 


9.47±0.49 


17.44±0.86 


Observed 


58 


54 


107 


219 



TABLE XI: Summary table by lepton flavor content, of background estimates, tt predictions and observed 
events in data corresponding to an integrated luminosity of 2.8 fb -1 for the >2 jet bin before the Ht and 
the opposite lepton charge requirement events. The uncertainties are the sums in quadrature of the statistical 
and systematic errors. The last column is the total dilepton sample obtained as the sum of the ee, and e/j, 
contributions. 



Nj ot > 2 tt Sample per Dilepton Flavor Category 



Source ee fi/j, e[i tl 

WW 3.54±0.63 3.65±0.65 7.50±1.28 14.70±2.47 

WZ 1.75±0.23 1.01±0.14 1.68±0.23 4.44±0.57 

ZZ 0.83±0.65 0.74±0.58 0.47±0.37 2.04±1.59 

14-7 0.62±0.41 O.OOiO.OO 1.45±0.58 2.07±0.78 

DY^ tt 2.97±0.75 3.29±0.84 6.68±1.67 12.94±3.22 

DY^ee + fifi 20.33±6.00 8.04±2.73 1.76±0.72 30.13±8.54 

W+jet fakes 9.66±3.00 18.67±5.77 54.67±16.59 83.00±22.90 

Total background 39.71±8.73 35.40±7.36 74.22±17.13 149.33±28.19 

tt (cr = 6.7 pb) 31.25zhl.52 32.69±1.59 74.62±3.58 138.56±6.61 

Observed 58 68 143 269 
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DATA (2.8 fb 1 ) 
(Entries : 684) 
Syst. Uncert. 
I tt (n = 6.7pb) 
i | Diboson 
H DY -> ee+p.p. 

□ DY -> TT 

■ Fakes 



*4o "sf Isf MO 




DATA (2.8 fb 1 ) 
(Entries : 342) 
Syst. Uncert. 
I Itf (a = 6.7 pb) 
I | Diboson 
H DY -> ee+tip. 

□ DY ^ TT 

■ Fakes 




DATA (2.8 fb 1 ) 
(Entries : 342) 
Syst. Uncert. 
I Itf (n = 6.7pb) 
I | Diboson 
H DY ee+trtr 

□ DY -> TT 

■ Fakes 



40 60 SO 100 120 140 160 180 200 
Dilepton invariant mass (GeV/c 2 ) 



DATA (2.8 fb 1 ) 
(Entries : 342) 
Syst. Uncert. 
I Itf (o=6.7pb) 
I | Diboson 
H DY -> ee+p.p. 

I DY -> TT 

■ Fakes 



60 SO 100 120 140 160 180 200 
Missing E T (GeV) 




200 300 400 
Hr(GeV) 



FIG. 5: Background and top quark signal predictions overlaid on the data for 0-jet events in 2.8 fb -1 . From 
top left to bottom right: Two leptons transverse energy spectrum, the dilepton invariant mass, the event Ifx 
and Ht- The hatched area represents the uncertainty in the total background estimate. 



TABLE XII: Summary table by lepton flavor content of background estimates, tt predictions and observed 
events in the final sample of events with > 2 jets passing all candidate selection criteria, for an integrated 
luminosity of 2.8 fb -1 . The uncertainties are the sums in quadrature of the statistical and systematic errors. 
The last column is the total dilepton sample obtained as the sum of the ee, /i/i and efj, contributions. 



tt Signal Events per Dilepton Fit 


ivor Category 




Source 


ee /i/i 


e/j, 


11 


WW 


2.16±0.38 2.42±0.42 


4.79±0.80 


9.37±1.51 


WZ 


0.94±0.15 0.68±0.11 


0.59±0.10 


2.22±0.33 


ZZ 


0.65±0.51 0.64±0.50 


0.23±0.18 


1.51±1.18 




0.23±0.25 O.OOiO.OO 


O.OOiO.OO 


0.23±0.25 


DY^ TT 


1.67±0.32 1.76±0.34 


3.87±0.72 


7.29±1.34 


DY— > ee + 


11.81±2.16 5.32±1.23 


1.36±0.60 


18.49±2.73 


M^+jet fakes 


3.91±1.28 9.34±3.05 


20.90±6.43 


34.15±9.51 


Total background 


21.37±3.14 20.16±3.64 


31.73±6.78 


73.26±11.30 


it (a = 6.7 pb) 


28.80±1.41 31.24±1.52 


70.15±3.36 


130.19±6.21 


Total SM expectation 50.17±4.25 51.40±5.00 101.88±10.06 203.45±17.33 


Observed 


41 55 


99 


195 
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-•- DATA (2.8 fb 1 ) 
(Entries : 438) 
'///. Syst. Uncert. 
I tt (<J=6.7pb) 
I I Diboson 
H DY -> ee+iiii 

I DY — > TT 

■ Fakes 



60 80 100 120 140 *^60 180 200 
Lepton P T (GeV/c) 



DATA (2.8 fb 1 ) 
(Entries : 21 9) 
Syst. Uncert. 
I 111 (n = 6.7pb) 
I I Diboson 
H DY -> ee+iiii 

I DY TT 

■ Fakes 





-•- DATA (2.8 fb 1 ) 
(Entries: 21 9) 
'///. Syst. Uncert. 
I Itf (tj=6.7pb) 
I | Diboson 
H DY -> ee+ii|i 

I DY TT 

■ Fakes 



20 40 60 80 100 120 140 160 180 200 
Dilepton invariant mass (GeV/c 2 ) 



> 
C5 



g 30 
UJ 

20 



Njet=l 



X Test : 0.451 



DATA (2.8 fb 1 ) 
(Entries : 21 9) 
'///. Syst. Uncert. 
I 111 (o=6.7pb) 
i I Diboson 
H DY -> ee+iiii 

I DY TT 

■ Fakes 



60 80 100 120 140 160 180 200 
Missing E T (GeV) 




FIG. 6: Background and top quark signal predictions overlaid on the data for 1-jet events in 2.8 fb -1 From 
top left to bottom right: Two leptons transverse energy spectrum, the dilepton invariant mass, the event Ifr 
and Ht- The hatched area represents the uncertainty in the total background estimate. 



precision measurements like the Z cross section 
measurement: 

<7 M -(iso) = 6.40 ± 0.75(stat) ± 0.49(syst) pb. 

The luminosity error is common to all of these 
subsamplcs and is not cxplicitcly quoted. 

X. CONCLUSIONS 

We present a measurement of the ti cross sec- 
tion at the Tevatron in a sample of data corre- 
sponding to an integrated luminosity of 2.8 fb -1 
collected by the CDF II detector. Using events 
with two leptons, large missing energy and two 
or more jets we select a sample with a signal 
top quark pairs almost a factor of two larger 
than the background. The contamination from 
SM sources is checked in lower jet multiplicity 
samples. From the excess of data over the back- 
ground predictions we measure: 

a ti = 6.27 ± 0.73(stat) ± 0.63(syst) ± 
0.39(lum) pb, 



or 

a ti = 6.27 ± 1.03(total) pb, 

consistent with the NLO standard model calcu- 
lation of 6.7^o'g pb. Yields in the ee, /i/i and e/j 
final states are in agreement with the predictions 
from lepton universality. 
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-•- DATA (2.8 fb 1 ) 
(Entries : 538) 
'///. Syst. Uncert. 
I tt (<J=6.7pb) 
I I Diboson 
H DY -> ee+^ti 

I DY — > TT 

■ Fakes 



20 40 60 80 100 120 140 160 180 200 
Lepton P T (GeV/c) 



DATA (2.8 fb 1 ) 
(Entries : 269) 
Syst. Uncert. 
I 111 (n = 6.7pb) 
I I Diboson 
H DY -> ee+^ 

I DY TT 

■ Fakes 





-•- DATA (2.8 fb 1 ) 
(Entries : 269) 
'///. Syst. Uncert. 
I Itf (n = 6.7pb) 
I | Diboson 
BDY^ ee+^ 

I DY -> TT 

■ Fakes 



20 40 60 80 100 120 140 160 180 200 
Dilepton invariant mass (GeV/c 2 ) 



DATA (2.8 fb 1 ) 
(Entries : 269) 
Syst. Uncert. 
I 111 (o=6.7pb) 
' I Diboson 
H DY -> ee+^ 

I DY -> TT 

■ Fakes 
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Missing E T (GeV) 




FIG. 7: Background and top quark signal predictions overlaid on the data for 2-jet events before the Ht and the 
opposite lepton charge requirement in 2.8 fb -1 . From top left to bottom right: Two leptons transverse energy 
spectrum, the dilepton invariant mass, the event Ifr and Ht- The hatched area represents the uncertainty in 
the total background estimate. 
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'///. Syst. Uncert. 


— 








□ tf (a = 6.7 pb) 










I I Diboson 
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DATA (2.8 fb 1 ) 
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I DY — > TT 

■ Fakes 
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FIG. 8: Background and top quark signal predictions overlaid on the data for top quark DIL candidate events 
in 2.8 fb -1 . From top left to bottom right: Two leptons transverse energy spectrum, the dilepton invariant 
mass, the event Ifr and Ht. The hatched area represents the uncertainty in the total background estimate. 
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TABLE XIII: Summary table by jet multiplicity bin of background estimates, tt predictions and observed 
events in data corresponding to an integrated luminosity of 2.8 fb _1 . The uncertainties are the sums in 
quadrature of the statistical and systematic error. The last column contains the candidate events with Ht > 
200 GeV and opposite sign lepton cuts applied. 



Control Sample and Signal Events per Jet Multiplicity 



Source 


jet 


1 jet 


> 2 jet 


H T + OS 


WW 


142.87±12.57 


40.01±4.24 


14.70±2.47 


9.37±1.51 


WZ 


11.65±0.74 


11.72±0.52 


4.44±0.57 


2.22±0.33 


ZZ 


8.79±6.78 


4.08±3.14 


2.04±1.59 


1.51±1.18 


Wy 


28.03±7.11 


8.14±2.27 


2.07±0.78 


0.23±0.25 


DY^ TT 


3.89±0.58 


17.87±2.99 


12.94±3.22 


7.29±1.34 


DY^ ee + fifi 


40.04±5.77 


40.83±7.28 


30.13±8.54 


18.49±2.73 


VK+jct fakes 


72.49±19.02 


94.23±26.16 


83.00±22.90 


34.15±9.51 


Total background 


307.76±35.84 


216.87±32.46 149.33±28.19 


73.26±11.30 


tt (a = 6.7 pb) 


0.67±0.06 


17.44±0.86 


138.56±6.61 


130.19±6.21 


Total SM expectation 308.43±35.87 234.31±33.28 287.89±34.70 203.45±17.33 


Observed 


342 


219 


269 


195 




FIG. 9: Dilepton candidate events (black point) by jet multiplicity. The stacked histogram represents the 
background contribution and the tt contribution for an assumed a t t = 6.7 pb. The hatched area is the 
uncertainty in the total background estimate. 



